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Abstract 

This paper studies the large time existence for the motion of closed hyper¬ 
surfaces in a radially symmetric potential. In physical, this surface can be 
considered as an electrically charged membrane with a constant charge per 
area in a radially symmetric potential. The evolution of such surface has 
been investigated by Schniirer and Smoczyk (Evolution of hypersurfaces in 
central force fields, J. Reine Angew. Math. 550 (2002), 77-95). To study 
its motion, we introduce a quasi-linear degenerate hyperbolic equation which 
describes the motion of the surfaces extrinsically. Our main results show 
that the large time existence of such Cauchy problem and the stability with 
respect to small initial data. When the radially symmetric potential function 
v = 1, the local existence and stability results have been obtained by Notz 
(Closed Hypersurfaces driven by mean curvature and inner pressure, Comm. 
Pure Appl. Math. 66(5) (2013), 790-819). The proof is based on a new 
Nash-Moser iteration scheme. 
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1. Introduction and Main Results 

Let £ be an oriented smooth closed manifold of dimension n , and (A f,g) 
be the Euclidean space, i.e. M. = R n+1 and g be the Euclidean metric. 
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Consider a smooth family of immersions F : [0,T] x E —> A4, we define an 
action integral of the form 


S(F)= [ T JC (F)-V(F)dt, 

Jo 

where /C is the kinetic energy and V = V + J is the total inner energy, V is 
a radially symmetric potential energy and J is the inner pressure. 

We fix a reference measure dp, on E with a smooth density function defin¬ 
ing a mass distribution on E. Then the total kinetic energy of all the points 
of the surface is 


K(F) = i j |9 t F| V 

We denote dp t as the induced surface measure of the induced metric 
g = F{t)*g at time t, ip denotes a smooth, radially symmetric function (re¬ 
flecting the presence of a central force) depending on s := The radially 
symmetric potential energy of the hypersurface F(E) is defined by 

V(F) = j v(s(F))dfi t , 


where v is defined by 


u(s) 



( 1 . 1 ) 


where 77 : R + —> R is a smooth function. 
The inner pressure is defined as 


J(F) 


~P log( 


Vol(F) 

Vol 0 
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where p > 0 denotes a parameter which determines strength of the influence 
of the inner pressure compared to the surface tension, Vol(F) denotes the 
enclosed volume of the surface .F(E). The initial enclosed volume VoIq as 
well as the constant p are included for scaling reasons. This inner pressure is 
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motivated by that of an ideal gas with constant temperature, i.e. proportional 
to VoF 1 {F). One can see [17], |18] for more details on the inner pressure. 
Then the action integral is 

S(F) = i f f \d t F\ 2 diidt- f f v(s(F))dntdt + p f log ( V ^p )dt. 

J 0 «/E J 0 «/E JO 0 

The Euler-Lagrange equations of functional S(F) is 

VaAF = f^F) ( -H(F ) + VF - (1.2) 

where H(F) denotes the mean curvature of F(E), v denotes the outer unit 
normal, |F| denotes the absolute value function of F, Vg t denotes the covari¬ 
ant derivative along F, i.e. 

VdAF = d tt F a + f a AF)d t F^d t F\ 


with f^ 7 being the Christoffel symbols of g. Here and in the sequel, we use 
the Einstein summation convention. 

Another setting of above problem is that closed hypersurfaces moves in 
Riemannian manifolds with density v(s). To our knowledge, Gromov jb] 
first studied the manifolds with densities as “nnn-spaces”, and mentioned 
the natural generalization of mean curvature in such spaces. Morgan and 
” Jjj considered the corresponding mean curvature H u = 

When the density in 


his collaborators 
H- < Vcu, v >, where the density is denoted by e' 


n tM 2 

the Euclidean space R n+1 is the Gaussian density (^-)^e Borisenko 

and Miquel [2| studied a flow of a hypersurface driven by its mean curvature 
associated to the Gaussian density. It is obviously that the Gaussian density 
is a special case of the density v(s) = exp(—| f* ^1dw ). 

The study of mean curvature flow can date back to Brakke [I], who in¬ 
troduced the motion of a submanifold moving by its mean curvature in arbi¬ 
trary codimension and constructed a generalized varifold solution for all time. 
Huisken j§, 9j showed that the mean curvature flow has much abundant and 
complicated behaviour, e.g. singularity and asymptotic behaviours. Schniirer 
and Smoczyk [2o[ studied the asymptotic behaviour of mean curvature flow 
in central force fields, where they chose the radially symmetric potential as 
(11.11) . For complex geometry, Chen and Tian jsj] used the mean curvature flow 
in a Kahler-Einstein surface to show the symplectic property being preserved 
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as long as the smooth solution exists. Chen and Li [4] produced holomorphic 
curves from a given initial symplectic surface by the mean curvature flow 
method. LeFloch and Smoczyk [l3| established a hyperbolic mean curvature 
flow, which is a strickly hyperbolic equation when the tangential part van¬ 
ishing. Then they obtained the local existence and singularity for such flow. 
Meanwhile, He and Kong [7j also introduced a hyperbolic mean curvature 
flow and established the local existence and nonlinear stability for this kind of 
flow. After that, a hyperbolic mean curvature flow similar to the well-known 
Hamilton’s Ricci flow was introduced by Kong and his collaborators jll[ 12]. 
and they obtained the short time existence and lifespan result. Recently, A 
new kind of more physical hyperbolic mean curvature flow was established by 
Notz |h|, its motivation is closed hypersurfaces moving driven by mean cur¬ 
vature and inner pressure. Meanwhile, the local existence of smooth solutions 
and the stability with respect to initial data was obtained. 

We notice that the dpi -term of equation (II.2ft prevents reparametrization 
of (1 1.2 [) (see j^] for the method of reparametrization) to remove the degener¬ 
acy of such a quasilinear equation. Hence a suitable approximation method 
should be explored. One of main results in this paper is the large time 
existence of equation (jl .2(1 . 


Theorem 1.1. Let F 0 : £ —>■ A4 = R n+1 be any smooth immersion with 
VoI(Fq) = VoIq > 0, and initial velocity F\ £ T(Ff TM). Assume that the 
radially symmetric potential function v has the form U.l\) . and 


0<|^|<M P , 0 < p < oo. (1.3) 

w 

Then there is e > 0 and a smooth family of immerisions F : [0, —> A i = 

_R n+1 solving the Cauchy problem \1.2 j) with small initial data F(0, •) = eF 0 
and d t F( 0, •) = eF\. Here T is a positive constant and T(FfTM) denotes 
the space of smooth sections in a vector bundle Fq TM, and e is a positive 
sufficient small parameter. 


We remark that we construct a small amplitude smooth solution, which 
can exist on time interval [0, ~^=] for a positive constant T. Here e is a 
sufficient small positive parameter measures the nonlinear effects. 

Next result gives the stability and uniqueness of solutions for equation 
(II.2p with respect to the initial data in the Euclidean space R n+1 . 
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Theorem 1.2. Let s > max{2, |}. Assume that equation ( Ej| ) with two 
different small initial data (FUFj) and (F 0 ,Fi) has two different solutions 
F and F in := {F e C? T : |||F||| S , T < R < 1}. If 

lll-Po — -foil |s, T + lll-Fl — A|||s,T < C(e), 


then we have 


\\\F-F\\\- S:T <0(e). 

In particularly, when F 0 = F 0 and Fj = F\, the solution of equation HI. A) 
on the time interval [0, is unique, i.e. F — F. Here O(e) denotes an 
infinitesimal and one can see the definition of the norm 111 • 11 | S) t in the section 

3. 


The organization of this paper is as follows. In Section 2, we derive some 
conservation laws on the solution of equation fll.2j) . In section 3, we estab¬ 
lish the existence of weakly linear hyperbolic system which will arise in the 
linearization of equation (11.21) . Section 5 is devoted to solving equation (jl.2[) 
by means of a new Nash-Moser iteration scheme (For general Nash-Moser 
implict function theorem, one can see [10, 15, Ebl |. In the section 5, we give 
the proof of the stability problem with respect to the initial data and the 
uniqueness of solution for equation (1 1.2 [) . In the last section, we discuss the 
extension of our main results to more general case, i.e. (A i,g) be an oriented 
smooth complete n + 1 dimensional Riemannian manifold. 
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2. Conservation Laws 

In this section, we derive three conservation laws, which are satisfied by 
the solution of (j 1.2 [) . More precisely, they are the total energy conservation, 
the momentum conservation and the interior momentum conservation. 
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Let X be a killing vector field on Xi = R n+1 . We define the total energy 
by 

^ j \d t F\ 2 dp- j v(s(F))dp t + p\og( V ^p ), 

1 / S »/ S 0 

the momentum with respect to J of a solution F of (11.2p by 


M x (F(t, ■)) = J < dtF,X(F) > dp, 
the interior momentum with respect to a vector field Y (diva^Y = 0) by 


Qy(F(V))= [ < d t F, F*Y > dp. 
J s 

For the total energy conservation, we have: 


Lemma 2.1. Let F : |0,T) x E —» Xi be a solution of \1.2\) . Then we 
have 

(1) . E(F(t, •)) = E(F(Q, •)) for allte[0,T). 

(2) . Let d t F T be the tangential part of dtF. Then we have 


d t (\\d t F\ 2 + v^) =< d t F, Vv > ^ + 


P 


dp 


dp Vol(F) 


< d t F. o > + v(divd t F T )^jL. 


dp 


dp 


Proof. The proof of (1) is simple, so we omit it. Now we derive the equality 
in (2). Using (ll.2jh we derive 




i T\dpt , _ ^ CZ T7" O'. dpt ! ^ „.dpt 


= < d t F , Wd t d t F > +v(divd t F )—— + v < d t F , Ho > —— + d t v 

dp dp dp 

= — v < d t F, Ho > -\-v< d t F, pF > < Q t p_ \/ v > 

dp dp dp 


P 


■ u > < d t F, o > ~~~ T v(divd t F T )^jd- + v < d t F, Ho > ^ + d t v ^ 
V ol(r) dp dp dp dp 


= < d t F, Wv > ^ + 


P . o t-j dpt /j- a oT\ dpt 
d,a ' l/crf(F) <d,F ' V> ~d)i + v ( dwd ‘ F 


where we use the form of u in the last inequality. 


□ 
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Lemma 2.2. Assume that there is a positive constant c such that 


n 

vdp t > cVol(F n + 1 ). 


Then there exists a positive constant K depending on Vol 0 , c, p and E 0 such 
that 

Vol 0 e~^ < Vol(F) < K, 


-\d t F\ 2 dp + [ vdp t < E 0 + plog(-^—). 

! ./v VoL 


Proof. This proof is simple, and can also be found in 17], so we omit it. □ 


For the momentum conservation, we have: 


Lemma 2.3. Let F : ]0,T) x £ —> Ai be a solution of lil.2\) . Then we 
have 

(1) . M x (F(t, •)) = M x (F(0, ■)) for all t e [0, T). 

(2) . Let d t F T be the tangential part of d t F. Then we have 


d t < d f F, X >= vdivX T 


^+<Xv^,X> 

dp dp 


P 


V ol{u) 


<u,X > —. 

dp 


Proof. The momentum conservation (1) can be obtained by integrating (2) 
with respect to dpo and dt and using that 



y a 9 \ 

< v,X > dp t = — 1 9=0 


d n 


Vol(F q ) = 0. 


So we only prove (2). 

We denote (j) q as the local flow of X, which is by definition an isometry, 
and set F q = (f> q o F. Then we have 


d t <d t F,X> 


< Vd^tF^X > + < d t F,S7d t X > 

-v <Hv,X> X > + < Vu^, X > 

dp dp dp 


P 


Vol(u) 


< v, X > 


dp t 

dp 
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T dnt <9, 1 _ fJ rrr ^d/j,t <9, , tT7l ^dfji t 


= vdivX , 

dp o q 

+ <Vv^-,X> 

dii 


g =0 log (dllt(Fq)) 
p 


d^ + d q 

_ y dpt 

Vol(u) ’ dii 


q=0 v{s{F q )) 


d^ 


= vdivX 


T^P't , „ VT dflit 


p 


dii^ <Wv d!i' X > + Vol(u) 


< v, X > 


dfk 

dii 


Furthermore, since A' is killing, we get 


— | g =o \og(dn t (F q )) = 0, 
°q 


— \ q=0 v(s(F q )) = 0. 


□ 


Lemma 2.4. Let F : [0,T)xE —> M. be a solution of H1.P\) and d t F T 
be the tangential part of d t F. Then we have 

(1). Q Y {F(t, •)) = Q y (F(0, •)) for all t e [0,T). 
d t < d t F,F*Y >= \div dlx (\d t F\ 2 Y). 


Proof. This proof is similar with Proposition 2.4 in 18J, so we omit it. □ 


3. Analysis of linearized weakly hyperbolic systems 


This section is to discuss the linearized system of (1 1.2 j) . which is a weakly 
hyperbolic linear systems. More precisely, it is obtained by decomposing with 
respect to time dependent subbundles into a system of coupled linear wave 
equations and linear ODEs. Here we prove the local existence of smooth 
solution of an abstract system for convenience. 

Firstly, we give some definitions of the norm and some inequalities. The 
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. Let T : Q —> E be a k- 
and F be the Frechet space 


following setting is inspired by the work of 
dimensional Riemannian vector bundle over E 
C°°([0,T] x E,f2) of smooth time dependent sections of Q. Assume that 
we have an atlas of coordinate charts ( x a , U a ) of £ such that a = 1, ■ • • , J, 
%a{U a ) = i?i(0), and the sets x“ 1 (i?i(0)) cover £, where Tu(0) C R n is a 
ball with radius 1. Meanwhile, for each such chart there are smooth time 
dependent local sections A = 1 ,•••,</, and r^ a \ d = l,-- - , d", of 
(d' + d" = k) defined on the domain of the chart that together form a 












basis of the fiber over each point in U a . For any other chart (xp, Up) with 
U a HUp ^ 0, we assume that the (p) and (p ), p e U a C\Up, span the 
same space and are bases for this space. The spaces spanned by the v^ and 
the t^ are orthogonal. For the specihc (fixed) coordinate chart, we omit 
the index (a) for convenience. Let he F. In each coordinate chart we take 
the decompose 

h = h± + hy = u a va + 

For any s € [l,oo], H s (F?i(0), R d ) xJ be the Sobolev space of the set of 
functions (u A a ^). The corresponding norm is 

J d' 

ni. = EE 11 ^(o) 11 H s (-Bi (0)) • 

o=l A=1 

For v, we use the norm 


J d! 

wi. = EE II ^o ° 3; a 1 ||H s (B 1 (0))- 

q=1 A= 1 

C'([0,r];H-(R 1 (0) ) R‘ , ') XJ ) denotes the function spaces with the norm 

l J d' 

ll‘ U llc'([0,T];H s (Bi(0),R d ') xJ ) = SU P EEE ||^M(^)||h s (Bi(0))- 

[°’ T ] i =o 0=1 A =1 

:= q =0 C z ([0, T\; H s_ *(Ri(0), R' ) xJ ) denotes function space with the 
spacetime norm 

2 J d! 

inilt = su p ii H —* ( j b 1(0)) - 

[°T] j =0 o=l A= 1 

We similarly define || ■ ||c s as the spatial C s -norm and the L 2 -Sobolev norm 
of order s in space and time. It is obviously that 

IM|h s (Bi(o)) < II m IIc 0 ([o,t];H s (Bi(o))) — IIMIUt- 

We make use of the following inequalities, which can be proven using 
classical methods of calculus; see, for example, the book of Chapter 13 in 
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[24| . Note that while versions of these estimates hold in all dimensions, as 
presented here the estimates are dependent on the dimension of £ being n. 
Generic constants are denoted by Co, ci,..their values may vary in the same 
formula or in the same line. 

Sobolev inequality: For s > | and u 6 H^R^O)) we have 
IM|l°°(.Bi(0)) < ||m||h s (Si(0))- 

Product inequality: For u, v G IP^R^O)) fl H s (i? 1 (0)) we have 

||ww||h s (Bi(0)) < ||w||h s (Bi(0))||^||l oo (Bi(0)) + |M|l°°(Bi(0))|MIh s (Bi(0))- 

In particular, if s > we have 

||?H|h s (Bi(0)) < |M|h s (.Bi(0))|M|h s (Bi(0))- 

We denote the spatially weighted Lebesgue spaces by L^(R n ), which is 
equipped with the norm 



Then the Fourier transform is an isomorphism between L^(R ra ) and H s (R n ). 
For T > 0, s > max{|, 2} and 0 < R < 1, we define 

B^ r := {ueC s T = n* =o C z ([0,T];H s - z (R 1 (0),R d, ) xJ ) : |||n||| s , T < R < 1}. 

We say that 

h = u a u a + r d r d 

is a solution of a weakly hyperbolic linear system if in each local coordinate 


chart (x a ,U a ) we have 

d tt u A - v(F)L a u a - v(F)N A r - v(F)Q A u A = h A , (3.1) 

d tt r d - v(F)M d r - v(F)P d u = h d , (3.2) 

with the initial data 

u A (0) = u A , d t u A ( 0) = u A , (3.3) 

r d (0) = r A , d t r d ( 0) = r A (3.4) 
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where v(F) = exp (—j /) 2 dw ) and 

L a u a = a Al i(t, x)didjU A + a Al (t, x, ua, x d )diU A + a A (t, x, ua , r d )u A ,(3-5) 
N A r = + n A (t,x,u A ,r d y + n A0 (t,x,u A: r d )d t r d 

+n A1 (t,x,u A ,r d )Y f b t^)j r (ff) d l- 1 -' ( 3 - 6 ) 

/j=i ■' s 

Q a u = q^(t,x,v A ,T d )diU B + qz(t,x,v A ,T d )u B + qz°(t,x,v A ,T d )d t u B 



M d r = m d {t,x,u A ,r d )r l + mf(t,x,v A ,r d )d t r\ (3.8) 

P d u = Ppj (t, x, u A , r d )djU B + pf{t, x, u A , T d )d t u B . (3.9) 

We assume all the coefficients in (13. II) - (13.21) and u A and r d to be smooth 
functions on x a (U a ), and supp tfp),, SU PP c fp) B c -/(5i(0)), a Aij (t,x) = 
a Ajl (t,x) and a Alj (t,x) satishes 

f y u e 2 <«< ,y ij e 2 . veeR", ( 3 . 10 ) 


for some positive constants p 0 < p\. On the other hand, we need that the 
operators are coordinate invariant under coordinate transformation on £ and 
under a change of basis between different {y^i T d^) an d pp). This 
means that the external force W = h A v i + h d T d is an element of Q. 

For the coefficient of the operators L A , N A , Q A , M d and P d , we define 
the local norm 

j J 

l L \s = ^2[L]s,a = YXY ||o Ay ||h s (Bi.(0)) + Y ll aA *llH s (Bi(0)) 

a=l o=l A,i,j A,i 

+ Y ll aj4 ||H s (Bi(0))), 

A 

J J 

[n]s = y^[^]s,o = Y(Y ii^iIh^b^o)) + Y : \\ n t Hh^bpo)) 

q.= 1 a=l A,i,j A,i 

+ Y H n rf°llH s (Bi(0)) + Y H nA1 |lH s (Bi(0)) + Y II^(o)i IIH s (Bi( 0)))j 

A,d A A,j 
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J J 

[Q\s = £[Q]s,a = £(£ IkgllH^BtCO)) + £ II?^bI|h s (B 1 (Q))) 

a =1 a= 1 A,B,i A,B 

+ £11^11^(0)) + £ |k A1 ||H s (Bi(0)) + £ 11 cf Q)j 11 H s (Bi(0))), 
A,B A A,j 

J J 

[M] s = £[M] SiQ = £(£ ||mf|| H s (Bi(0)) + £ ||^°||h s (Bi(0))) ; 
a=l a=l d,i d,i 


[^]s - £[-P]s,a - £(£ \\Pb IIh s (B!( 0)) + £ 11 11 H s (Bn. (0))) - 

a=l a=l B,d,j B,d 

Assume that 

[L]s + [dtL] s + [A] s + [Q]s + [M]s + [P]s < c 0 , (3-11) 

det(< v A , v B >) > ci, det(< r dl r t >) > a. (3.12) 

We give some estimates on v = anc ^ Td = ®dF. The proof of the 

following results is similar with the results in page 26, 49 of 17j], so we omit 
it. 

Lemma 3.1. Let F € B S RT . Assume that ||^||co + UrUco < R for some 
R > 0 and Ii3.12\) holds. Then there exists a constant c\ such that 

11 1/45 11 fr (£ 1 ( 0 )) < ci(l + |Ml.ff s (£i(o))), 

11^ 11 IP (£ 1 ( 0 )) < Ci(l + 11'T11 (£i( 0 ))) , 

|||^ AS |||h 3 ([0,T]x£i(0)) < Ci(l + ||^||if 3 ([0,T]x£i(0))), 

11 111 iP([0,T] x£i(0)) < Ci(l + ||r||iP([0,T]x£i(0))), 

Ill'll. < Ci(l + ||F|| s+ i), IIMIU < Ci(l + |||F||| s+ i), 


ll^ll. < ci(l + £ II^FlUi), f° r i > !, 

1=0 

ll zyj4S |lc°(£ 1 (0)) < c i, \\T dl \\ cO(£i( 0 )) < Ci- 

From above result, we can see 

IMU + IM S+l + II^HIs + IM| S +1 + IIM.+1 < C 2 • 

The following result is the main result in this section, which states the 
existence of linear system 03.ip - fj3.2p . 


12 








Proposition 3.1. Assume that \3.1(h) - \3.1&) holds. Letho,hi G C°°(E,n) 
and F G B S RT be given. Then system A3. 1\) - FFE) has a unique smooth solu¬ 
tion h on [0, T\ x £ with /i(0) = ho and dth(0) = hi for some given external 
force W = h A ui + h d r d G C°°([0, T] x £, 12). 

Before giving the proof of above result, we carry out some priori estimates 
on the solution ( u A , r d ) of system fj3.ip ~ fj3.2j) in local coordinates. We remark 
that the following energy estimates does not depend on the compact property 
of the spatial domain. We also do not use the integral on spatial variable by 
part. The main idea of proof of Lemma 3.2-3.3 divides into two steps. The 
first step is to find a suitable differential inequality with respect to the time 
variables. The last step is to integrate the time variable by part. Meanwhile, 
the integral on spatial variable is taken. 

Lemma 3.2. Let (u A ,r d ) be a smooth solution to system A3.1\) - A3Ad\) with 
initial data A3.3\) - \3.4\) . F G B S RT is given. Assume that A3.10l) - K3.12\) holds. 
Then there holds 

[ T [ [A - c 3 v(F ) - c 4 n^P < d,F + d t F, F > v(F)}e~ xt 

Jo TBi(0) T I 

x [( u t) 2 + a Alj diU A djU A + ( u A ) 2 + ( r d ) 2 + ( r d ) 2 ]dxdt 

< f e~ M ((uf) 2 + a Mj d i u£d j u£+ (u£) 2 + {rf) 2 + (r$) 2 )dx 
Asi(o) 

+c 5 f f e~ xt [{h A ) 2 + ( h d ) 2 + {d ir d ) 2 ]dxdt. (3.13) 

Jo Jb^O) 

Proof. Taking the inner product of the linear system (13.11) with 2 e~ xt uf and 
using (13.5p . we have 

d t [e~ xt {{uf) 2 + (F)a Azj diU A djU A )\ + A e~ xt ((uf) 2 + v(F)a Atj diU A djU A ) 
-2 e- xt d i (v(F)a Aij ufd j u A ) 

= - 2v{F)e~ xt {d i a Aij )u?d i u A 

-2 (d t v(F) + diV^e-^a^d^dj^ 

+2 v(F)e- xt a Ai ufd i u A + 2 v(F)e~ xt a A ufu A 

+2 v{F)e~ xt ufN A r + 2 e~ xt v{F)ufQ A u A + 2 e~ xt u A h A . (3.14) 

Taking the inner product of the linear system (13.21) with 2e~ xt r d , we have 
dt[e~ Xt (rt) 2 ] + \ e ~ xt {r d t ) 2 
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Since we have 


2e xt v(F)r d M d r + 2e xt v(F)r d P d u + 2e xt r d h d . (3.15) 


d t [e~ xt {u A ) 2 ] + A e~ xt {u A ) 2 = 2 e~ xt u A u A , (3.16) 


d t [e~ xt ((r d ) 2 + (d ir d ) 2 )} + A e~ xt [{r d ) 2 + (<%r d ) 2 ] = 2 e~ x \r d r d + (3.17) 


and 

MF) = <: SF,F > «(F), 

summing up 03.141) - (13. 171) . it holds 


dt[e~ xt ((u?) 2 + v(F)a Aij d i u A d j u A + ( u A ) 2 + (rf) 2 + ( r d ) 2 )] 
+Ae- At [( M f) 2 + v{F)a Aij d i u A d j u A + (u A ) 2 + (rf) 2 + (r d ) 2 } 

-2 e~ xt d i {v(F)a Aij u^d j u A ) 

= v(F)e- At (a t a^)a i M A a i M A - 2u(F) e - At (d i a^K A d i « A 

< diF + d t F, F > v{F)e~ xt a Aij d i u A d j u A 

\F\ 

+2 v(F)e- xt a Ai ufd iU A + 2 v(F)e~ xt a A ufu A 
+2v(F)e~ xt u A N A r + 2 e~ xt v(F)ufQ A u A + 2 e~ xt u A h A 
+2 e~ xt v(F)r d M d r + 2 e~ xt v(F)r d P d u + 2 e~ xt r d h d 
+2 e~ x \u A ut + r d r d + d^rf). (3.18) 

Note that v{F) > 1. Using Cauchy inequality, by (I3.6I) - (13.9|) . Q3.10[) - (l3.1ip . 
(13.18P and Lemma 3.1, we derive 


d t [e~ x \(ui) 2 + v(F)a Aij d i u A d j u A + (u A ) 2 + ( r d ) 2 + (r d ) 2 )] 
+A e~ xt {{u?) 2 + v(F)a Aij d i u A d j u A + {u A ) 2 + (rf) 2 + (r d ) 2 } 
-2 e- xt d i {v(F)a Aij u A d j u A ) 


< c 3 v(F)e- xt [(u?) 2 + ( diU A ) 2 + (a ,-« A ) 2 + (u^) 2 + (rf) 2 + (r d ) 2 ] 

.^(l-^l 2 ) , a r a r fx cn„-At|7a „.A\2 , „.A\2i 


+c 4 n- 


l*T 


< diF + F > v(F)e~ xt [{diU A Y + (d jU A y] 


+4 e- x \d ir d t y + c 6 e- At [(h A ) 2 + (h d ) 2 }. 


d \2 


-\t[/uA\2 


d\ 2l 


(3.19) 
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Then inequality (j3. 19j) leads to 


d t [e~ xt {{uf) 2 + a Aij d iU A d jU A + ( u A ) 2 + (rf) 2 + {r d ) 2 )} 

+ [A - c 3 v{F) - + d t F, F > v{F)]e~ xt 

\ F \ 

x[{u?) 2 + a A ^d iU A djU A + {u A ) 2 + {rf) 2 + {r d ) 2 } 

-2 e~ xt d i {v{F)a Aij ufd j u A ) 

< 4e~ xt {d ir f) 2 + c 6 e~ xt [{h A ) 2 + ( h d ) 2 ]. 

Thus integrating both side of above inequality on [0,T] x £>i(0), we obtain 


/ e~ xt [{uf{T )) 2 + a Aij d iU A {T)d jU A {T) + {u A {T)) 2 

Jb^o) 

+{rf{T)) 2 + {r d {T)) 2 }dx 

+ [ [ [A - c 3 u(F) - c A n^p- < d t F + dtF, F > v{F)}e~ xt 

Jo Jb 1 ( 0 ) \ F \ 

x[{u?) 2 + a Aij d iU A d jU A + {u A ) 2 + {rf) 2 + {r d ) 2 }dxdt 


< 


'Bi(O) 


e~ xt [{uf) 2 + a Ai W iU A d jU A + (■ u A ) 2 + {rf) 2 + {r d 0 ) 2 }dx 


+4 



e- xt {d ir f) 2 + c 6 


'o Jb i(0) 

This completes the proof. 



e xt [{h A ) 2 + {h d ) 2 ]dxdt. 


o Jb i(o) 


□ 


Lemma 3.3. Let {u A ,r d ) be a smooth solution to system H3.1\) - (33^) with 
initial data \ 3. A) - /)3.4\ ). F 6 B S R T is given. Assume that i \1.3\) and 13.1CK) - 
1)3. 1 A) hold. Then there exists a positive constant A > C 7 such that 


(A - c 7 ) 
< c 8 [ 



e~ xt [{uf) 2 + {0 iU A ) 2 + {u A ) 2 + {rf) 2 + {r d ) 2 ]dxdt 


,A\2 


A\2 1 f„d\2 


0 Jb 1(0) 

e~ M [{u A ) 2 + (diU A ) 2 + {uf) 2 + {diuf) 2 

d\2 1 2 1 /Q „d\2] 


Jb 1(0) 

+(r d ) 2 + {d t r d 0 ) 2 + {rf) 2 + {d ir f) 2 ]dx 
+c 9 [ [ e~ xt {\h A \ 2 + | dih A \ 2 + \h d \ 2 + Ojh 2 )dxdt. 


(3.20) 


0 JBi( 0 ) 
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Proof, we introduce auxiliary functions u A and r d , which satisfy 


d tt u A — d t u A — u A = h A : (3.21) 

u A { 0, x) = u A , d t u A ( 0, x) = u A , 

and 

d tt r d - d t r d - f d = h d , (3.22) 

f d ( 0,x) = Tq, 9 4 r d (0,x) = r d . 

Let 

u A = u A - u A , (3.23) 

f d = r d _ f d , (3.24) 

then it follows from f!3.If) - (13.211 and (I3.2ip ~ fl3.22p that 

d tt u A - v(F)L a u a - v(F)N A f - v(F)Q a u a = g u (3.25) 
■u A (0,a;) = 0, d t u A (0, x) = 0, 


d tt f d - v{F)M d f - v(F)P d u = g 2 , (3.26) 

r d (0,x)=0, d t r d (0, x) = 0, 

where L A , N A , Q A , M d and are defined in fj3.5p ~ fj3.9jl . and 

g x = v(F)L a u a + v(F)N A f + u(F)Q a m a - d t u A - u A , (3.27) 
g 2 = v{F)M d f + v(F)P d u - d t f d - f d . (3.28) 

Using the similar deriving process with (13.1911 . from (I3.25jl - fl3.26p we get 

d t [e~ xt {{u A f + v(F)a Aij d i u A d j u A + ( u A ) 2 + (rf) 2 + (f d ) 2 )] 
+\e~ xt [(u? ) 2 + V (F)a^'c^c^ + (h^) 2 + (rf) 2 + ( f d ) 2 ] 

—2 e~ xt di(v(F)a Alj ufdjU A ) 

= v(F)e~ xt (d t a Aij )d i u A d j u A - 2v{F) e - x \d l a Mj )ufd i u A 
+ n V ^P < diF + d t F, F > v(F)e~ xt a Aij d i u A d j u A 
+2 v(F)e~ xt a Al ufdiU A + 2 v(F)e~ xt a A ufu A 
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+2v(F)e~ xt ufN A r + 2 e~ xt v(F)ufQ A u A + 2 e~ xt v{F)r d M d r 
+2 e~ xt v{F)r?P d u + 2 e~ xt (u A u? + r d r d + 8#%$) 

+2e~ M u A gi + 2e~ xt r d g 2 . (3.29) 

In order to avoid an extra loss of derivatives of integrating (j3.29|) . we need to 
rewrite (I3.29j) . By (13.27|) ~ (13.28j) . the last two terms in (I3.29P can be rewritten 
as 

e xt {ufgi + r d g 2 ) 

= d i (v(F)e~ xt a Alj d i u A d t u A ) - d t (v(F)e~ xt a Aij d i u A d j u A ) 
-v{F)e- xt d i a Aij d i u A d t u A - Xv(F)e~ xt a Aij d i u A d j u A 
+v(F)e~ xt d t a Aij d i u A d j u A + v(F)e~ xt a Aij diU A d t djU A 

< diF - d t F, F > v(F)e~ xt a AlJ diU A djU A + v(F)e~ xt a Al diU A u t 

+v(F)e~ xt a A u A Ut + e~ xt u?(v{F)N A r + v(F)Q A u A - d t u A - u A ) 
+e~ xt r d (v(F)M d f + v(F)P d u - f d - f d ). (3.30) 

Note that v(F) > 1. Inserting (13.301) into (13.29f) . integrating in [0, T] x i?i(0), 
then using Cauchy inequality we obtain 

[ e~ xt [(uf{T)) 2 + v(F)a Aij d i u A (T)d j u A (T) + (u A (T)) 2 + (ff (T)) 2 

Jb x { o) 

+(r d (T)) 2 ]dx + [ [ [A-| ^P-<d i F + d t F,F>\v(F)] 

Jo Jb!( o) l-H 

xe~ xt [(u?) 2 + v(F)a Aij d i u A d j u A + {u A f + ( r d ) 2 + ( r d ) 2 ]dxdt 
< [ e~ xt v(F)(d i u A {T)) 2 dx 

+c w n f f e ~ A *l ^11^1 ^ < diF ~ 9 t.F- F > | v(F)(diU A ) 2 dxdt 

Jo isi(o) \F r 

+c u ( [ e~ xt v(F){(diuf) 2 + (diU A ) 2 + (uf) 2 + (u A ) 2 dxdt 

Jo J Bi(0) 

+ (rf) 2 + (f d ) 2 + (d.rf^iM. (3.31) 

So it follows from (I3.23p - (l3.24p and (I3.3ip that 

ft l*~\ !! TW- <d , F + d t F,F>\v(F)] 

Jo Jb x { 0 ) \F\ 
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xe~ xt [(uf) 2 + v{F)a Aij d i u A d j u A + ( u A ) 2 + (rf ) 2 + ( r d ) 2 }dxdt 

< [ e~ xt v(F)(d i u A (T)) 2 dx+ [ [ [. \-\?!QS)- < d i F + d t F,F>\v(F )] 

J -Bi(O) Jo J -Bi(O) F I 

xe - At [(« t A ) 2 + u(F)a A ^ A 6 tu A + (w A ) 2 + (rf ) 2 + (f d ) 2 ]cM 

+c 10 n [ f e~ xt \< d t F - d 4 F, F > \ v(F)(d l u A ) 2 dxdt 
Jo Jb!(0) f I 

+c n [ T [ e- xt v{F)[{diU A ) 2 + (dpi A ) 2 + (uf ) 2 + (u A ) 2 

io J -Bi(O) 

+(rf ) 2 + (f d ) 2 + (diff) 2 ]dxdt. (3.32) 

Multiplying (13.21ft and (13.22ft both side by 2 e~ xt u A and 2 e~ xt r d , respectively, 
then integrating on [0, T] x -Bi(O) and summing up two inequalities, we have 

[ e~ xt [(u A (T )) 2 + ( d t u A (T )) 2 + (f d (T )) 2 + (air d (T)) 2 ]dx 

Jb^o) 

+X [ [ e~ xt [(u A ) 2 + (d t u A ) 2 + {r d ) 2 + (d t f d ) 2 ]dxdt 

Jo Jb^o) 

< cu [ e ~ xt [(u A ) 2 + «) 2 + (rf ) 2 + (rtfjdx 

Jb^o) 

+C 13 f j e~ xt (\h A \ 2 + | h d \ 2 )dxdt. (3.33) 

Jo ifii(O) 

Differentiating (I3.21[) with respect to x, by the similar process of getting 
(137331) . we derive 


IB i(0) 


+A 


e- At [(a i h A (T)) 2 + ( dAu A (T )) 2 + (a,f d (T)) 2 + {dtdif d (T)) 2 ]dx 



e-^Kd ^) 2 + (dAu A ) 2 + ( d t r d ) 2 + (d^f d ) 2 ]M; 


- A\2 


d\2 




0 ^Bi(O) 


< c 14 / e- At [«) 2 + «) 2 + (a,<) 2 + (a,<) 2 ]dx 

J Bi(0) 


+cis f e~ xt [(r d ) 2 + (rf) 2 + (d.rf) 2 + (ftrf) 2 ]^ 

J Bi (0) 

+c 16 [ [ e~ xt {\h A \ 2 + |a,h A | 2 + |/r d | 2 + | dih d \ 2 )dxdt. (3.34) 

Jo isi(0) 
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On the other hand, we notice that 


F £ Br jT , ||F|| l ~ < \\F\\n“ and H-FUcodo^jH 8 ) < 111-^11 It, s f or s > 1. 
So by fl 1.31) . we have 

( ^| 2 


< 


d i F + d t F,F> | < |||F|||^ + 1 (||F|| c o ([ 0iT];HS+1) + ||F|| c 1 ([ 0 jT]; h S) ) 

< 2IIIFI 


|f + ?T < ^ P+2 < 1, 


v(F) < e%F\°\ Pd °- < e 2 <P+V(p+i) 


2(p+l) 

s,T 


_n_ 

< g2(P+ 2 )(p+l) 


ij2(p+l) 


(3.35) 


Hence by fl3.32j) . there exists a A such that 


_n_ 

A > e 2(p+2) (p+n 


i?2(p+l) 


> 1, 


and 


_ n _r 2 ( p +i) 

(A - e 2(p+2) (p+i) ) 


"_R2(p+1) 


'0 JB i( 0 ) 


e Xt [{uf) 2 + e 2(p+2) (p+!) a Alj diU A djU A 


+{u A ) 2 + (rf) 2 + (r d ) 2 ]<M 


d\2 , ✓ d\2l 


7 , -_R2(P+1) 

< g2(P+ 2 ) (p+1) 


e xt (diU A (T)) 2 dxdt 


JBi{ 0) 

__»2(p+l) 

+(A - e 2(p+2) (P+D ; 


'0 ^Bi(O) 


_n_ o2(p+l) 

e _A *[('uf ) 2 + e 2 (p+ 2 , (p+i)' K a^di^djii A 


+(u A Y + (rf ) 2 + (f d ) 2 }dxdt 

. , _2_p2(p+l) I 

+Ci 7 ne _At e 2(p+2) (p+i) / 


(diU A ) 2 dxdt 


+Cae ^krrX w '[ [ 

Jo </Bi(0) 

+ (ff ) 2 + (f d ) 2 + (^ff) 2 ]dxdt. 


-A\2 

'0 ^Bi(O) 

.-At 17 0 ~A \2 


e- Ai [(^wf ) 2 + (a iU A ) 2 + (tif ) 2 + ( u A ) 


~A \2 


(3.36) 


On the other hand, inequalities fl3.33p ~ fl3.34p give the control of four terms 
in the right hand side of f!3.36p . Thus substituting (I3.33p ~ fl3.34p in (I3.36p . we 
obtain (I3.20p . This completes the proof. □ 
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In what follows, we plan to obtain the estimate for 11 u A \ | | Sj r and 11 \r d \ | | Si t 
by considering the equation of time space derivatives of u A and r d . For conve¬ 
nience, we recall the norm |||w||| Si t = sup [0iT] E?=o ELi Ea=i II^^IIh-Wo))- 
For any multi-index 7 6 Z 2 + with |y| = s, applying dJ to both sides of (13. 1 j) 
and (13.2j) to get 

d tt &lu A - v(F)L A 8 1 i u A - v{F)N A 3 1 i r - v(F)Q A 8 1 i u A = h A , (3.37) 
d tt &lr d - v(F)M d d]r - v(F)P d & y i u = h d , (3.38) 

where the nonlinear terms 


h A = 


d?h A + d 81 (v(F)a Aij )d i d j d 82 


u 


Sl+S2=7' 

, |si|>l 



E 

d?{v{F)a M )didl*u A 


ai+S2=7> |si|>1 




E 

d 8l (v{F)a A )D S2 u A 

+ e wn 

nf)did 82 r j 

ai+S2=7i |si|>l 


Sl+S2=7, |si|>1 


E 

d 8l (v(F)n A )D S2 r j 

+ E af'(y(F)nf)^rf 

Sl+S2=7> bl|>l 

d- l (v(F)„ A1 )J2 

31+52=7, |si|>1 


E 



ai+S2=7> |si|>l 

P=1 

E 


E 

d’ A (v(F)q$)d l d!*u B + E &‘'(v(F)q A B )dfu B 

Sl+S2=7! |si|>1 


si+s 2 =7, |si|>1 


E 

d-‘(v(F) g f)d»u B 


si+S2=7> |si|>1 

d!'(v(F)q A1 ) E I 



E 

^ d ?( C {P)Bufl3)) d V 

(3.39) 

Sl+S2=7> |si|>1 

P=1 ^ 

E 



and 

q = d?h“+ V d!‘(v(F)mi)d-V+ J2 Sl'(v(F)mf>)dtH 

S1+S2=7, |si|>1 si+S 2=7, |ai|>l 

+ E 9" l (v(F)P$))a j dpu B 

«1+S2=7, |si|>1 
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(3.40) 


+ V d?(v(F)Pg)d?v?. 

si+s 2 =7, |si|>1 

For convenience, we denote all spacial derivatives of u A and r d of the order 
s by a column vector of m(s) components and n(s) components 

u 7 s = (d^u A , dl~ 1 d 2 u A , ..., d s 2 u A ), 
r J = (, d s 1 r d ,d s 1 - 1 d 2 r d ,...,d%r d ). 

Putting (13.3711 - 03.381) together corresponding to all 7 with j'y] = s, we have 

d tt u s - v(F)L a u s - v(F)N A r s - v(F)Q A u s = h A , (3.41) 
0 tt r s - v(F)M d r s - v{F)P d u s = h d , (3.42) 

where 

L A u s = a Al ^didjU s + a Al diU s + a A u s , 



M d r s = m d r s + fh do d t r s , 

P d u s = PgdjUs + P£°d t u s . 


with all the coefficients in L A , N A , Q A , M d and P d are (m x m)-matrices 
and h A and h d are m- vectors given by 

U A — (h A h A h A \ T 

^7 = (^fs,0,...,0)) rfs- 1,1,...,0); ■ ■ • ) ^(0,0,...,s)) • 

Lemma 3.4. Let ( u A ,r d ) 6e a smooth solution to H3.1\) - T3T^) with initial 
data \3.3\) - \3.f\) . F E B S RT is given. Assume that HI. 5 1) and H3.1(h) - T3.1'A) 
hold. Then for s > 2, there holds 

11 M |\ S ,T + || M I |s,T < Ci 9 (| |rt 0 | | S +1,T + 11^11 |s+l,T + I | r o| |s+l,T + I | r l||s+l,T 

+lll'i1ll.,T +111*111., t ), (3.43) 

where <79 depends on T. 


21 















Proof. The proof is based on the induction. For s — 1, Lemma 3.3 gives the 
case by choosing a suitable A. We assume that (13.4711 holds for all 1 < j < s, 
and we prove that (I3.47P holds for s + 1. Since (I3.4ip - (l3.42p has the same 
structure with (13. 1 p ~ (13.2p . the result in Lemma 3.3 can be used directly. Note 
that the vector D s u A and D s r d satisfy (I3.4ip - fl3.42p . By (I3.20p . we have 


(A — c 2 o) 



e- xt \(dlut) 2 + (9* +1 u A ) 2 + (dtu A ) 


,5+1 A\2 


45 A\2 


'0 JB i(0) 

2 1 /Qs„d\2i 


+(d s i r d y + ( dyy]dxdt 

< c 2 i 

J B\ (0) 

3S+1 2 1 ins„d \2 


1 e-»mu A y + (<rv) 2 + (<9> A ) 2 + (<rv) 2 

■JB i(0) 

+ (<9*r 0 d ) 2 + (dt +l r d 0 ) 2 + (, dtr d ) 2 + (<9* +1 r?) 2 ]dx 

T 


+<* E / / 

|^,| =s ^0 9Bi(0) 


e- At (|<| 2 + + |/*T + |«| 2 )M.(3.44) 


„A|2 


By (1373511 . dT39H - (l3T(l . we drive 


12 


VI + + \h d \ + 

< \dth A \ 2 + | W 2 + |aVl 2 + \dt +l h d \ 2 

+c 23 (| W| 2 + \d s iU A \ + 19?" VI + 19* r d | + larVfl). (3.45) 

So choosing a A > c 2 o + c 2 4 +1 large enough, it follows from (I3.44p and (I3.45P 
that 


(A — C 20 — C 24 ) 

'0 J -Bi(O) 
+{d s i r d ) 2 + (d s ir d ) 2 ]dxdt 



e- M mu?Y + (ar l u A Y + ( 9 > 


A\2 


< C 21 


'Si(O) 


e- x, mu$Y + (9' +1 « Y? + + (S‘ +1 «f) : 


+Wh) 2 + (cf + V 0 J ) 2 + (afrfi 2 + (af +1 r?) 2 ]& 


e-^da'A- 4 ! 2 + lay 1 /!- 4 ! 2 + layy + layvi 2 )**, 



+C25 

' 0 4BP0) 

which implies that 


II^V||l2( B i(0)) + II^ +1uA |Il 2 (Bi( 0)) + ll^^ r<i ||L 2 (Bi(0)) + II^VlIlA^O)) 
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< c 26(||Wo IIh 3+1 (Bi(0)) + II U 1 1 |Ih s + 1 (Si(0)) 

+ ll r ollH s + 1 ( J Bi(0)) + ll^l IIh s + 1 (Si( 0)) + H^llH'+^BitO)) + 11 h d 11 h s + 1 (_Bi (o))) • (3.46) 

Furthermore, we can anolv did? +] ~ 3 to both sides of (I3.1I)-(I3.2D for 2 <j< 
s + 1, then deriving a similar estimate with (13.461) . We conclude that 

s+l 

^(ll^t^ +1_JMA |lL a (B 1 (0)) + ||^^ +1 "V|| La ( fll (o))) 

3=0 

< C27{\\Uq IIh s + 2 (Bi(0)) + II u ^IIh s + 2 (Bi( 0)) + II r 0 IIh s+2 (Bi(0)) 

+ ll r l IIh s + 2 (Bi(0)) + II^‘ 4 |Ih s + 2 (Bi( 0)) + ||h d ||H»+2(Bi(0))). (3.47) 

Summing up the estimate (13.471) over all coordinate charts, and we notice 
that \\u\\ s is equivalent to ||w|| s + ||5^u|| s , hence (13.431) can be derived by 
(13.471) . This completes the proof. □ 

Note that h = u A v A + rd Td and W = h A u A + h d Td. So by Lemma 3.5, 
direct computation gives the following result. 

Lemma 3.5. Let h be a smooth solution to weakly linear hyperbolic system 
( 1 3. with initial data h( 0) = h 0 and d t h(0) — h\. F G B S RT is given. 
Assume that M.3\) and Ii3. ICh) - \3.lZfy hold. Then for s > max{2,|}, there 
holds 


lll^llls, T < C 28 (||ho|| s+ i )T + ||^l||s+l,T + |||W|||s,t)- (3.48) 

Proof of Proposition 3.1. The proof is based on a standard fixed 
point iteration. Let h = u A u A + r d Td and W = h A v A + h d Td■ We consider the 
following approximation system of f l3.1l) - fl3.2D as 

dttuf m+ 1 ) - v(F)L A uf m+1] - v(F)N A r {m) - v(F)Q A uf m) = h A , 
dur d [m+ 1 ) - v(F)M d r {m) - v{F)P d u [m) = h d , 

with initial data 

w m+i(°) = u t r i+ 1(°) = r t 

d t u A +1 (0 ) = u A =< hx, i/ B {0) > u AB (0) - u° < dtv c {G), z'b(O) > u AB (0 ) 
-Tq < d t T d (0), z'b(O) > u AB ( 0), 

d t r d m+ M = r d =<h 1 ,Ti(0)>T ld (0)-u A <d t v A (0),T l (0)>T lk (0) 
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-r 0 m <a t r m (0),r i (0)>^(0). 


Above system is a linear wave equation coupled with an ODEs. By a standard 
fixed point iteration and similar estimates in Lemma 3.5, we can prove that 
it has a local smooth solution on [0,T]. For the existence of solution for 


general linear wave equation, one can see [171 . 1181 22, [23J for more details 


4. Proof of Theorem 1.1 

In this section, we construct a solution F(t, x) on [0, of equation (11.2[) 
by a Nash-Moser iteration scheme. Rescaling in (j 1.2 j) amplitude and time as 

F(t, x) eF(y/et, x), e > 0, 


we obtain the following system 


d ti F = s 


d/j, 


v(sF)(-H{sF) + scp(sF)F 


-£V(eF)V(^) + 


P 




2 ' ' Vol(eF) J ’ 

with initial data 

F(Q,-) = F 0 , d t F(Q,-) = F 1 . 

Introduce an auxiliary function 

F = F — F 0 — t,F \, 

then the initial value problem (14.1[) - (I4.2[) is equivalent to 

d tt F = s^^vieiF + Fo + tF^i-HieiF + Fo+tF,)) 
d(i 

+s(F + F 0 + tF 1 )ip(s(F + F 0 + tFi)) 

-e\(e(F + F„ + tF,))V( l?+ F °„+ tF ' ? ) 

+-=-^- )u, 

Vol{s{F + Fo+tFjy 


(4.1) 


(4.2) 


(4.3) 


(4.4) 


with zero initial data 


f( 0,0 = 0, a t F(o,-) = o. 


(4.5) 
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We define 


A f(F) := d tt F-£ n - 1< ^v(£(F + F 0 + tF 1 ))(-H{£(F + F 0 + tF 1 )) 

a/i 

+e{F + F 0 + tF 1 )tp(e(F + F 0 + tF,)) 

+ + tFO)v( | 7 + Y tFl12 ) 

+ ^o((e(F+F c , + iF 1 ))*‘ / ’ ^' 6) 

then our target is to construct a F(y/et, x ) E B^ T on the time interval [0, T] 
such that A/"(F) = 0. Here E [0,F], Thus we obtain the existence of 
smooth solution for (j4.4[) with initial data fl4.5jl . In fact, we treat the initial 
value problem (14.4D - (14.51) iteratively as a small perturbation of the initial 
value problem for a weakly linear hyperbolic system. Linearizing nonlinear 
equation (14. 4 j) . we obtain the linearized operator 


dpJ\f(F)h := d tt h — £ n 1 ^-u(e(F + F 0 + tiq))[(A < h, v > +\h\ 2 < h, v > 

UjLXj 

— < ViL, h >) + £hip(£(F + Fq + tF,) + e{F + Fq + tFi)(dp(p)h) 


-£ 2 (dp(p)hV( 

P 


I F + Fq + tF\\~ 


) - £ 2 ip(£(F + F 0 + fFi))V < F + F 0 + LFi, h > 


< h, v > d)i t ]v 


Vol{e(F + F 0 + tFi)) 2 
-e n -^ v (e(F + F 0 + tFi))[—H(e(F + F 0 + tF\)) 
+e(F + F 0 + tFi)ip(£(F + F 0 + tFi)) 

-eM4F + f„ + (F,)) V (!L±4±tLJi) 


P 


](divh T + < h,u > H)p 


Vol(£(F + F 0 + tFi)) 

+£ n ^ V (£(F + F 0 + tF 1 ))[-H(e(F + F 0 + tF,)) 

an 

+e(F + F 0 + tFi)ip(e(F + F 0 + tFi)) - £ 2 cp(£(F + F 0 + LF\))V( 
P 


\F + F 0 + LFi| 


Vol(£(F + F 0 + fFi)) 


]<^(e(F + F 0 + fFi)) < |F + F 0 + fFi|, h > v 
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“ £n_1 ^( £ ( F + Fo + tF ^~ H « F + F o + tF 1 )) 
+e(F + F 0 + tFi)(^(e(F + F 0 + LFi)) 


-£ 2 if(s(F + F 0 + tFi))V( 


I F + Fq + tF\ | 


P 


](V < h, v > -e 2 < h, diF > h dl d d F ), 


Vol(e(F + F 0 + tFi)) 


(4.7) 


where h l3 denote the second fundamental form, dp denotes the Frechet deriva¬ 
tive to F and we use the following relations 

-dpH = A <h,u> +\h\ 2 <h,v>-< WH , h T >, 
dpdpLt = (clivh 1 + < h, v > H)dfjL t , 
dpu = V < h, v > - < h, d{F > h ld d d F, 

dpV^siyF + Fq + tFi)) = — £ip(s)v(£^F + Fq + tFi)) < \F + Fq + tF \|, h > . 


For the nonlinear term, by (j4.4jl and (14.71b direct computations show that 
77(h) := A f(F + h)~ Af(F) - dpAf(F)h. (4.8) 


Since the exact form of nonlinear term (14. 8 p is very complicated, here we 
does not write it down. In fact, we notice that the solution of (14.4(1 is to be 
constructed in B^ T and F G B^ T , so we have 


ll^llls,T — lll^llls,T’ V 9 >2. 


Thus we only need to obtain the lowest exponent of h for convenience. Then 
by the product inequality, we can obtain the following result, i.e. Lemma 
4.1. Note that 


v(s) = exp(—^ / } -^-dw), d^t = Jdei{g~j), 
z 11 w v 


9iFx...x d n F 
\d\F x...xd n F 


Vol(F) 


V oIq + 



< d t F, v > d[i t ds. 


Direct computation shows that the lowest exponent of F is n + 1 in (14. 4 p . 
which is determined by the lowest exponent of F in d/i t - This lowest exponent 
of h in (14.81) equals to the lowest exponent of F in (14.4ft . i.e. n + 1. 
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Lemma 4.1. Let F, h e B S RT . For any s > max{2, there holds 

IIIWWIIUj-^IIWIKt. («) 

where c 2g depends on the constant R. 

In order to prove the existence of smooth solution for the linearized equa¬ 
tion, we can follow jl8| to decompose h and the external force W (t, x ) into 
the normal and tangential parts. Then direct computation shows that the 
linearized equation has the same form as the weakly linear hyperbolic system 
(I3.1l) - (l3.2|h So by Proposition 3.1 and Lemma 3.5, we obtain the following 
result. 

Lemma 4.2. Let F e B S RT . For any s > max{2, |}, there exists a 
function h e B S RT such that 

dpAf (F) h = W(t, x), 

h{ 0,-)=0, d t h( 0,-) = 0. 

Moreover, there holds 

ll|ft|||.,T < c 3 „|||»'|||., T . (4.10) 

Next we introduce the smooth truncation function, one can see (2l| for 
more details. Let D Iq e C°°(R) such that = 0 for 9 < 0 and II^ —> / for 
6 —y oo. we introduce a family of smooth functions S(9') with S(9') = 0 for 
9' < 0 and S(9') = 1 for 9' > 1. For W G C^, we define 

U WHxl W(t,x) = S'Gar'l - \x\)W(t,x), x,x' G R n . 

For l — 0,1, 2,...,, by setting 

\x'\ := iV; = 2\ (4.11) 

then by the Fourier transform is an isomorphism between L^ and H*, it is 
directly to check that 

||n M _ N fF|| H n < c illSa ^ 1 - 8a ||^|| H * a , V Sl > s 2 > 0, (4.12) 

WB-\x'\-\x\W - WWns, <c SUS 2 N^- S 2 \\W\\ U s 2 , V 0 < < s 2 . 
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For convenience, we denote IIjv,-|x| by ILv r We approximate system (14. 6 1) . 
and get the following approximation system 


G Nl {F) := a«F- £ "- 1 n iV ^n(£(F + F 0 + tF 1 ))[-i/( e (F + Fo + tF 1 )) 
+e(F + F 0 + tF 1 )ip(e(F + F 0 + tFi)) 

-s*tp(e(F + F 0 + tF 1 ))V( |F + F ° + tFl12 ) 


+ 


P 


Vol(e(F + F 0 + tFi)) 


(4.13) 


where m = 0,1,.. 

The following Lemma is to construct the u l t h" step approximation solu¬ 
tion. 

Lemma 4.3. There exists a solution h l+1 of the initial value problem 

E l + dpM(F)h l+1 = 0, 
h l+1 (0, x) = 0, d t h l+1 (0, x) = 0, 

where E l satisfies 

E l = G Nl (F l ). 

Moreover, for s > max{2, |}, it holds 

|||^ +1 |IUt < c 3 i|||^||| Sj t. 


(4.14) 


Proof. Assume that a suitable ”0tli step” approximation solution of (14.13li 
has been chosen, which is W° 0. The ‘7th step” approximation solution is 
denoted by 


i =0 


Dehne 
E l = 


_n— Itt dpt ~fT7i l 


d tt F u - £ n ~ 1 U Nl ^-v(e(F l + F 0 + tFi))[—H(s(F l + F 0 + tFj) 
+e(F 1 + F 0 + tFi)ip(e(F l + F 0 + tFfi) 






(4.15) 


i -=/ 




■e 2 v(e(F + F 0 + tF,))V( 
P 


I F +F 0 + tF 1 \- 




W 


Vol(e(F +F 0 + tF',)) 




Then we plan to find the “7th step” approximation solution F . By (14.13|) . 
we have 

g(F l + h l+1 ) 

= d tt F l + duh l+1 - e^U Nl+ i^v(e(F l + h l+1 + F 0 + tF,)) 


x[-H{s(F l + h l+1 + F 0 + tF 1 )) 




+e(F + F 0 + tF,)(p(e 
-e 2 <f(e(F l + F 0 + tF,))\7( 

+___ P _] 

Vol(e{F l + h l + l + F 0 + tF,)) 
= E l + dpAf (F)h l+1 + R(h l+1 ), 


+ tF,)) 


| F+h l+1 + F 0 + tF,\' 


where 


(4.16) 


R(h l+1 ) = -e n - l Il Nl+1 [Af{F l + h l+1 ) - U Nl Af{F l ) - fiL,A f{F l )h l+A j . (4.17) 

By Lemma 4.2, there exists a solution h l+1 of 

E l + dpAf (F)h l+l = 0, 
h l+1 (0, ■) = 0, d t h l+1 ( 0, -) = 0. 

A similar estimate with (14.10p is derived as 

ll|ft i+1 |ll.,T<c 31 |||i3 , |||», 7 .. 

Furthermore, one can know from (14.161) and (j4.17[) that 

E l+1 = R(h l+1 ). (4.18) 

□ 


29 














For max{2, |} < s < s 0 < s, set 

sr-=s+ s -^, (4.19) 

s — s 

a l+1 := si - s i+ 1 = (4.20) 


By (14.191) - (14. 2 Op . it follows that 


s 0 > Si > ... > si > s i+ 1 >- 

Theorem 4.1. Equation ( |^.1[ ) with initial data (4-2) has a solution 

F = Fqo + Fq + Fit, (4-21) 

where F^ has the form 


F 00 = J2 h ^ C r• 


i=0 


Proof. The proof is based on the induction. For any l = 0,1, 2,..., we claim 
that there exists a constant 0 < d < 1 such that 


|||^ m ||U +1 ,T<d 2i <1, 

ll|£ m | 


U+i — ^ 

T‘ +1 e b^ t . 


-=0 


, 2 i+i 


We choose a hxed sufficient small F > 0 such that 


|F°| 


iso 


<1, lll£° 


|SQ 


< 1. 


For the case l = 0, by (I4.14p . we have 


I Sl 


,T < C3i|||£' 0 ||| S1) T < C 31 I | |-F°| | |s 0 ,T < 1- 


(4.22) 

(4.23) 

(4.24) 

(4.25) 

(4.26) 


It follows from (I4.9j) . (14.12!) . ,17ft and (14.18ft that 

lll^llkr < ||kk)llls 1 ,T 

< C32 ^- i ivk i iiih i iii: i k 


< c 3 3(2it|||B 0 ||| w ,t) 


n+1 


(4.27) 


30 




















It is obviously to see that fl4.2bp - fl4.27p gives fl4.22lh _n — (I4.23p ?-n by choosing 
suitable small e > 0. So we get F e B'^ T . 

Assume that (j4.22n ~ fl4.24p holds for 1 < i < l, i.e. 


III^ 1 IL,t < i, 

infill Si <d 2 \ (4.28) 

r g b£ t . 

Now we prove that (j4.22jl - fj4.24j) holds for l + 1. From (I4.14|) and (I4.28p . we 
have 


|||^ +1 ||U+i,T — C3l|||^||| Si ,T < C 31 f/ 2 < 1. 

It follows from (14.9p . (14.12j) . (14.14p . (14 .1 7p and (14.18p that 


pl~\~ 1 I 


\Sl + l 


,T < |||fl(ft' +1 )|||. I+1 ,T 


< c 3t e n -‘N. 


n— 1 /\ 7 -n+l 

l-\~ 1 


< c 3t E n -'N;^\ 


l~\~l I I I 71 -(-1 

11 ls i+1 ,T 

j^l 11 \n+l 


I si,T 


< c 34 


I E 


,«-ii I ibH- 1 ) 2 
I bi_i,T 


< 


< C3 5 (16£|||S 0 ||| So> t) 


(n+l) 


1+1 


(4.29) 


(4.30) 


We can choose a fixed sufficient small £ > 0 such that 

0< 16£|||£°||| S0 , r < 1. 

Thus we conclude that (I4.22p - fl4.23p holds. Note that F l = ^ =0 /d. So 
(M gives «. 

Therefore, we derive 

lim |||^|||;,t = 0 , 

l —1-00 

which implies that system (14.4p with zero initial data has a solution 

OO 

Foo = eC * T - 

i =0 

At last, by f!4.3p we obtain the solution of system (14. ip with initial data f!4.2p 
has a solution 


F — F oo + Fq + F\t. 


□ 
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5. Proof of Theorem 1.2 


This section is to study the stability problem of (1 1.2 in B^ r . Assume 
that there exist two solutions F and F of (14.11) with corresponding to two 
different initial data ( Fq, F±) and (Fq, Fi). Let F = F — F. Then we have 


d tt F = s' 


t -i dfit(F) - dnt(F) 

d/x 


v(eF) 


F I 


x -H(eF) + £ip{eF)F - £ 2 ip(sF)V{ L ^-) + 




Vol{£F) J 


+£ n - ld ^^-{v(£F) - v(£p)) 

an 


x \ -H(£F) + £<p(£F)F - eV(eF)V(^-) + 

+ £ n- 1 ^l v (eF) 

an 

f _ _ \F\ 2 

x -H(eF) + £y(£F)F - eV(eF)V( J — L ) + 


VolpF) 


P 


Vol(£F) 


u(F) 


u(F) 


(u(F) - u(F)) 


+ £ n—1 d Pi( F ) v( £ F)(-H{£F) + H{£F) + £(tp{£F)F - <p(,sF)F) 


dn 


1*1 




e 2 (v(eF)V(^-) - )) + 


P 


P 


VoI{eF) VoI(eF) 


HF),(5.1) 


with initial data 


Fo-F 0 - F (h Fi — F l — Fi. (5.2) 

Introduce an auxiliary function 

F* = F-F 0 - tFl , (5.3) 

then the initial value problem (15.If) - (I5.2[) is equivalent to a zero initial value 
problem. For convenience, we denote the right hand side of (15. ip by Af(F, F, F). 
Then equation (15.ip can be rewritten as 

d tt F* - J\f(F*, F,F) = 0. (5.4) 

We consider the approximation system of (15.4p as 

Q'(F*) := d tt F* - n JVi W(0, F, F). (5.5) 
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Then using the similar computation process with (14.161) . we have 


C e (F*) + E'(t, x ) + R(F*) - E\t , x) = 0, (5.6) 

where E'(t,x) is a function which does not depends on F*, 

C £ (F*) = d tt F * - n Nl d F .Af{0, F. F)F*, 

R(F*) = £ n ~ 1 U Nl (A f{F*, F, F) - Af(F, F, F) - d F *Af(0, F, F)F* (5.7) 

Direct computation shows that C e (F*) is a weak hyperbolic system, which 
has the same structure as (ID. ID - (13.21) . So by Proposition 3.2, there exists a 
solution F* of 


C £ (F*) + E'(t,x) = 0, 
Fq = 0, d t F* = 0. 

A similar estimate with (14.141) is derived as 


111 F* j j | s } T < C 3 g| | \E'\ | \s,T- 

Since the lowest exponent of F* in (I5.7j) is n + 1, so we have the following 
result. 

Lemma 5.1. Let F, F e B S RT . For any s > max{2, there holds 

II|W(F-)|||., t <c 37 |||F*|||J« t , (5.8) 

where C 37 depends on the constant R. 

Then by (15. 6 p and (15.8f) . we have 


IF* 


l,T < C 36 | | l-E'l | |s/,T < C 36 |||i2(F*)| || Sii (T 


I /-I* 11 \n +1 
T 11 ls;,T 


< c 38 £ n_ 1 A7 l+1 | 


< c 39 £ n - 1 ^ n+ 1 A'/"+ 1)2 | 

< ... 


IF* 


< C4o(8£|||F-|||„,t) < "+ 1 >‘, 
which combines with (15.31) . we have 


I (n+1) 2 
I Sl-2,T 


IFI 


I s u T < 


I-Pol I U,T + T\ 11 F x | || SI) T + C 40 (8£| 11 F* 11 \s 0 ,t) 


(n+l) ! 
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Thus we obtain 


|||F-^||| S!l T<|||Fo-F 0 ||| Si)T + T|||F 1 -F 1 ||| S/iT + C4o(8£|||F*||| so , T )^ ) 

which implies the stability result for equation (j 1 . 2 j) . 

In particularly, if two solutions F and F of (14.1 p has the same initial data, 
i.e. Fq = Fq and F\ — F\. Then we choose a suitable small e such that 

0 < 8e|||F*||| a0) T < 1. 

Then we have 

lirn \\\F- F||| ShT = 0, 

L - >-00 

which gives the uniqueness of solution for equation ( 1 1 . 2 j) with the initial data 
F 0 and iq. This completes the proof of Theorem 1.2. 


6. Further Discussion 

Our main results can be extended to more general case, i.e. (M,g) be 
an oriented smooth complete n + 1 dimensional Riemannian manifold. This 
section will discuss this case shortly. By Nash’s embedding theorem we may 
embed (M , g ) isometrically into an Euclidean space R A by f : A4 —> R A 
for some large k. Then we use the similar method of [lfl, 0] to derive an 
extrinsic form of the Euler-Lagrange equation (11.21) . 

We denote T m as the closest point projection to /(Af), it is defined on 
a neighborhood 

M = {x + y\xe f(M ), y e {T x f(M)) L , \y\ < 5(x)}, 

of f(M) and is smooth there. Here 5(x) denotes a positive smooth function 
on f(Ai). Since the second fundamental form h a p of Jv[ is normal to /(At) 
and / = T M o /, we get 


dadpf - TZpd.f = D A D B ^ M dpf A d a f B , (6.1) 

where D A and D B denote the derivative in the direction of the canonical 
basis vector e A and e B in R fc , respectively. 
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Let F : [0, T\ x £ —► M and F : [0, T] x £ —> f(M) with F = f(F). 
Note that / is an isometric embedding. Then by fl6.ll) . direct computation 
shows that 


n g F - ^ v (F)—f—u(F)-D A D B ^ M (F)(d t F A d t F B -^v(F)g^d i F A d j F B ) 
dg vol(r) dg 

- D 1 ^v(F) V (F)(r 1 oF~V(}-^-)) V (F) 

= D,(y e ,a t F - ^ v(-H(F ) + cpF- v>V(/L) + 

where 


□ g F = «9«F - ^ v (F)g i ^d i d j F - T^d d F). 


Hence F is a solution of (II.2ft if and only if F is the solution of 


-D A D B V M (F)(d t F A d t F B - ^j^ V (F)g ij diF A djF B ) = 0. (6.2) 

Next we follow the method of [jl| to extend equation (16. 2 p for F : [0, T] x 
£ —> M C R a which do not necessarily map to /(Ad). Let 'L^x(.F) be the 
projection onto the normal space of Q t = F(t, £). We notice that if F is close 
enough in C 1 to a family of immersions that map to /(Ad), then o F is 
also a family of immersions and u(^f m ° F) and v {^>m ° F ) is defined. Thus 
we will deal with the following equation 


u g F- d Jl~ v ^JL_ 
dv vol(F) 


V 


*n f i .D A D B ^ M {F){d t F A d t F B - ^vjidiF^F 3 ) 
t dg 

<P fii LB / ^^ ( r 1 oF- v (lh!!))P=0, (6.3) 


where v = o F), v = \l/ a xz/(\l /^ o F) and Ip = d^x^d/vi o F). 

The following result shows that if F maps to /(Ad) initially and the initial 
velocity is tangent to /(Ad), then F maps to /(Ad) for all time. The proof 


is the same with Lemma 4.3 in 18 , so we omit it 
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Lemma 6.1. Assume that equation H6.3\) with initial data F( 0,x) G 
/(A4) and d t F(0,x ) G Tf(o,z)/( M) has a smooth solution F on [0, T] x E. 
Then F(t,x ) G /(A4) and / -1 o F solves equation HI. 2 1) with initial data 
f- 1 o F( 0) and Df~ 1 (d t F(0)) on [0, T]xE. 

In what follows, one can use the same process of dealing with the Eu¬ 
clidean case to prove the existence of solution for equation (16.3p . Here we 
omit the details of proof. 
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